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Abstract 

We study the universal thermodynamic properties of systems consisting of many coupled oscilla- 
tors operating in the vicinity of a homogeneous oscillating instability. In the thermodynamic limit, 
the Hopf bifurcation is a dynamic critical point far from equilibrium described by a statistical field 
theory. We perform a perturbative renormalization group study, and show that at the critical 
point a generic relation between correlation and response functions appears. At the same time the 
fluctuation-dissipation relation is strongly violated. 

PACS numbers: 64.60.Ht,05.45.Xt,05.10.Cc 
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The collective behavior of many interacting elements generally leads to transitions and 
critical points in the large-scale and long-time properties of complex systems. This is well 
understood in the study of systems at thermodynamic equilibrium iL^^lsf. Non-equilibrium 
critical behaviors have been studied in a number of systems jj, Joj but remain a serious 
challenge. An important example for criticality far from thermodynamic equilibrium is the 
behavior of coupled oscillators in the vicinity of a continuous homogeneous oscillatory insta- 
bility or supercritical Hopf bifurcation. Such instabilities are important in many physical, 
chemical and biological systems 0, 0] . 

In this paper, we apply concepts of the theory of dynamic critical points to study the 
generic properties of systems of coupled oscillators in the thermodynamic limit. In partic- 
ular, we discuss linear response and two-point correlation functions defined for the oscilla- 
tor ensemble. Since the system is far from a thermodynamic equilibrium, the fluctuation- 
dissipation (FD) relation between correlation and response functions in equilibrium systems 
is broken. We show that a Hopf bifurcation represents a non-equilibrium critical point and 
study the universal behaviors characterizing its approach from the non-oscillating state. We 
apply field theoretic renormalization group (RG) methods and develop an RG procedure 
which is appropriate for the case of a spontaneously oscillating system. This RG is per- 
formed in an oscillating reference frame with a scale dependent oscillation frequency. The 
RG fixed points characterize the universal critical properties of locally coupled oscillators. 
We find that at the critical point of a Hopf bifurcation an FD relation is formally satisfied if 
the system is described within the oscillating reference frame. In terms of physical variables, 
the FD relation is strongly broken but a relationship between correlation and response func- 
tions appears. Even though our calculations are performed in a d = A — e dimensional space, 
we suggest that the main features of our results apply to Hopf bifurcations in general. 

The generic behavior of an oscillator in the vicinity of a supercritical Hopf bifurcation can 
be described by a dynamic equation for a complex variable Z which characterizes the phase 
and amplitude of the oscillations Q]. This variable can be chosen such that its real part is, 
to linear order, related to a physical observable, e.g. the displacement X{t) generated by 
a mechanical oscillator: X{t) = Re(Z(t))+nonlinear terms. In the presence of a periodic 
stimulus force F{t) = Fe~^'^^ with a frequency uj close to the oscillation frequency at the 
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bifurcation uq, the generic dynamics obeys 

dtZ = -{r + iuJo)Z -{u + tUa)\Z\^ Z + A-^e'^F{t) . (1) 

For F = and r > 0, the static state Z = is stable. The system undergoes a Hopf 
bifurcation at r = and exhibits spontaneous oscillations for r < 0. The nonlinear term 
characterized by the coefficients u and Ua stabilizes the oscillation amplitude for u > 0. The 
external stimulus appears linearly in this equation and couples in general with a phase shift 
6. In the case of a mechanical oscillator, the coefficient A has units of a friction. From 
the point of view of statistical physics, the Hopf bifurcation is a critical point and Eq. 
characterizes the corresponding mean field theory. Indeed, at r = and u = Uq and in 
terms of the amplitude X of the limit cycle Z{t) = Xe"*"^*, the system exhibits a power-law 
response |X| ~ where 5 = 3 is a mean field critical exponent. For frequency differences 

Icj — ciJol ^ A~^/'^|Fp/^|u + mal^''^, the response becomes linear with |X| ~ A~-'-|F|/|co' — tuol- 
In the presence of fluctuations, the critical point of an individual oscillator is concealed in 
the same way as finite size effects destroy a phase transition in equilibrium thermodynamics. 
However, a true critical point can exist in a thermodynamic limit where many oscillators, 
distributed on a lattice in a (i-dimensional space, are coupled by nearest-neighbor interac- 
tions. The combined system undergoes a dynamic phase transition at which all oscillators 
synchronize and an order parameter, which characterizes the global phase and amplitude of 
oscillators, becomes nonzero. Since at the critical point the correlation length diverges and 
oscillators become synchronized over large distances, a discrete model can on large scales 
be described by a continuum field theory which characterizes the universal features of the 



critical point 
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In the case of localh^ 



the complex Ginzburg-Landau equation 



coupled oscillators, this field theory is given by 



111 ] with fluctuations 



dtZ = -{r + iuJo)Z + {c + iCa)AZ 

-{u + iUa)\Z\^ Z + A-^e'^F + r] . (2) 

Here, Z(x, t) becomes a complex field defined at positions x in a (i-dimensional space and 
A denotes the Laplace operator. The coefficients c and characterize the local coupling of 
oscillators and the effects of fluctuations are described via a complex random forcing term 
r7(x, t). For a vanishing external field F(x, t) and in the absence of fluctuations, Eq. (0) is 
invariant with respect to phase changes of the oscillations Z —>■ Ze^'^. As far as long time 



and long wavelength properties are concerned, t] can be chosen Gaussian and white with the 
correlations (r7(x, t)?7(x', t')) = and (r7(x, t)?7*(x', t')) = AD6'^{'x — x.')6{t — t') which respect 
phase invariance. 

The linear response function Xaf3 and the two-point auto-correlation function Cap of 
this field theory are defined by (^„(x,t)) = / d'^x' dt' Xafsi^ - x',t - t')Fp{x.' ,t') + 0{\F\^) 
and Cq,/3(x — x',t — t') = {ipa{^,t)ipf^{x.' ,t'))c- Here we have expressed Z = ipi + iip2 and 
F = Fi + iF2 by their real and imaginary parts and (...)c denotes a connected correlation 
function. Because of phase invariance, these functions obey symmetry relations, e.g. Cu = 
C22, C21 = —Ci2- In the following, we focus for simplicity on the elements C = Cu and 
X = Xn which characterize the correlation and response of the observable X. 

It is convenient to eliminate the frequency ujq from Eq. Q by a time-dependent variable 
transformation Y = e*'^o*Z, H = e*'^o*A"^e*^F and ( = e*'^''*?]. This leads to the amplitude 
equation 

dtY = -rY + {c + iCa)AY-{u + iUa)\YfY + H + ( , (3) 

where the noise ( has the same correlators as rj. For the particular case Ca = and Ua = 0, 
Eq. (j2)) becomes identical to the model A dynamics of a real Ginzburg-Landau field theory 
with an 0(2) symmetry of the order parameter j^. The critical behavior of this theory at 
thermodynamic equilibrium has been extensively studied This leads, in this particular 
case, to a formal analogy between an equilibrium phase transition and a Hopf bifurcation. 
The correlation and response functions Cajj and Xai3 can here be obtained from those of the 
equilibrium field theory by using the time-dependent variable transformation given above. 
Since the theory at thermodynamic equilibrium obeys an FD relation, a generic relation 
between the correlation and response functions Cap and Xap appears. This special case 
provides a further example of an equilibrium universality class found in a non-equilibrium 

n 

dynamics with non-conserved order parameter It is the case e.g. for the model A dynam- 
ics of the real Ginzburg-Landau theory with a Z o sy mmetry [IS], even when the symmetry 
is broken by the non-equilibrium perturbations [l^, and for some of its generalizations to 
the 0{n) symmetry 15]. 

This raises the question to know whether the equilibrium universality class also charac- 
terizes the general case where Ua and are finite. Dimensional analysis reveals that for 
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d > 4 mean field theory applies. In this case, 



-ie 



+ 



2A [R - i{uj -no) R- i{uj + ^o) 
D D 



+ 



(4) 



where R = r + cq^, Qq = ^^o + Caq^ and where q and u are wave vector and angular frequency, 
respectively. 

For d < 4 mean field theory breaks down. We apply RG methods using an e expansion 



near the upper critical dimension (c? = 4 — e) Defining two real fields 0^ by F = (pi+i(p2 
Eq. (jS)) reads 



-Rap4>t3 - f4/3</>/307</>7 + + (a , (5) 
(r - cA)5af3 - CaASap and U^fS = Ud^p + Ua^^a/J, with €21 



where H = Hi + iH2, Rq/s 

— ei2 = 1 and e^j = for i = j. We introduce the Martin-Siggia-Rose response field 
and apply the Janssen-De Dominicis formalism to write a generating functional with 
action 



d'^X dt |D0Q0Q, — (pa [dt4>a + Ra 

-Uaf34>a(pf3(p-y(t>-y^ ■ (6) 

Using a Callan-Symanzik RG scheme , we define the renormalized theory such that 
its effective action is of the form ©. This requires to introduce a phase shift 66 and a 
frequency shift 6ujq between the bare fields 

Q^pi-Sujot) Zl^'^Z^ 0^(x, t) 

nap{-se-6ujot)zy'z^M^,t) . (7) 



t>a,<Pa) renormalized fields 



Here we have introduced Z-factors for the renormalization of the fields and the time 
(t" = Z^^t), and Qai3{6) denotes the rotation matrix by an angle 6 in two dimensions. 
We furthermore introduce dimensionless coupling constants g and Qa and a scale factor /i by 
(47r)-^/2^ and Un. = /i'^(47r) '^^'^ga- The bare and renormalized quantities are now re- 
lated depending on fx [2^. We define the correlation and response functions Gaf3 = C 
and 7q/3 = {4>a4>i3)c- They are related to the physical observables Cap and Xap via |2f 
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Capiy^^t^) = VLaa{,-UJQt) Z^ZlGap{yi,t) 



Xaf3 



;x,t°) = A-'naAo-uJot)iz^z.Y/'zl^,pi^,t) . 



(8) 



The dependence of the renormahzed parameters g, Qa and Ca on /i defines three /^-functions. 
Writing g = {g,ga,Ca), we have P {g,e) = fi{d^g)Q, where /3 = {P,Pa,Pc) and (9^)o denotes 
differentiation with fixed u^, and c°. 

To one-loop order in perturbation theory (see Fig. ^ for examples of Feynman diagrams 
of the theory), only Uq, r, g and ga are renormalized. The two non-trivial /^-functions are 
given by: 



p -eg -V 
Pa ^ -ega - 



a 



:{g^ - gl + '^ggaCa) -6gga 



(9) 



where V = 4i5/(47r)^. The RG fixed point corresponds to the values g* of the parameters g 
for which the three /^-functions are simultaneously zero. Since /3c = to one-loop order, one 
condition is lacking to fully determine g*. Choosing Cq as a parameter, we obtain g* ^ e/5V 
and gl ^ Cqc/SP. In order to determine completely the fixed point g*, we need to go 
to two-loop order in perturbation theory. To this order all parameters (apart from c and 
D), the fields and the time are renormalized and explicit expressions for all the Wilson's 
functions and Z-factors of the theory can be obtained. 

Only one fixed point exists which describes the universality class of Hopf bifurcations. It 
obeys c* = 0, (?* = and is infrared- stable. This fixed point is formally equivalent to the 
one of the real Ginzburg-Landau theory with 0(2) symmetry. As a consequence, we find 
z/ ~ 1/2 + e/10, 1] ^ and 2; ~ 2 + (61n(4/3) — 1) /50 which are the corresponding 

equilibrium critical exponents. Here, u denotes the exponent characterizing the divergence of 
the correlation length ^, z is the dynamic exponent and t] denotes the exponent characterizing 
the field renormalization 

The RG fiow in the vicinity of the fixed point however is defined here in a larger parameter 
space as the one corresponding to the 0(2) dynamic model. Furthermore, the effective theory 
discussed here is expressed in an oscillating reference frame with scale-dependent frequency 
and phase. Therefore the correlation and response functions Cap and Xa/s differ from those 
of the equilibrium model and additional universal exponents appear. We can derive generic 
expressions for these functions using the RG flow of all parameters in the vicinity of the 
critical point and employing a matching procedure [2^. For example, we flnd for ^ 1 
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and for stimulation at the effective oscillation frequency Uq^: 

pie{q) 



1 



(10) 



g2-^ 2Aeff 

where we have introduced the functions 9{q) ~ O^s + ctcsq'^^ + Pesl'^'^ and 7(g) ~ 7eff?'^^ 
and non-universal effective quantities denoted by the index "eff" . The universal exponents 
Ui ~ e/5 and uj2 — here are characteristic for a Hopf bifurcation. Similarly, we find 

expressions for the correlation function: 

. (11) 

and for the frequency dependence of the homogeneous mode g = in the regime (a; — 

Because of the formal analogy of the RG fixed point discussed here with the one of an 
equilibrium field theory, an FD relation appears exactly at the critical point and relates the 
functions Gai3 and •ya/s- Since the physical correlation and response functions Cap and Xai3 
can be determined from G^fs and using the scale dependent variable transformations of 
Eq. (jS)), a relation between correlation and response functions appears: 

cos6lcfrXii - sin6'efrXi2 = oa (^'^n + iL^fCi2) 

costei2 + sinteii = OA n {^fCu + iuCu) ■ (12) 

Here, Xa/s = x'a^ + has been separated in its real and imaginary parts. The relation 
(fT^ is asymptotically satisfied in the long time and wave-length limits at the critical point. 

The physical correlation and response functions Ca/s and Xa/3 however do not obey the 
equilibrium FD relation. The degree of this violation can be characterized by a frequency- 
dependent effective temperature Tgfj: 

y;ff(^) ^ ^ Cn{uj,q = 0) . . 

T 2kBT x'liiuj,q = 0) ' ^ 

Here, fc^ denotes the Boltzmann constant and T is the temperature of the system. We find 

that Teff/T ~ (cu — uJq^)^'^ diverges at the critical point with a universal exponent a. For 

the particular case = and Ua = 0, a = 1, while otherwise cr^ 1 — e/5 to first order in e. 

The power-law divergence of T^g as a function of frequency reveals a violent breaking of the 

FD relation. This divergence at the oscillation frequency has been experimentally observed 



on a single active oscillating system 



2l| 



We have shown that the critical point in d = A — e dimensions of locally coupled oscillators 
is formally related to the equilibrium phase transition in the XY model. From this analogy 
it follows that on the oscillating side of the Hopf bifurcation, the system of coupled oscillators 
exhibits long range phase order and coherent oscillations for d > 2. We can speculate how 
our results are modified in lower dimensions d. In analogy with the equilibrium XY model, 
we expect the phase order of the oscillations to vanish for d < 2, and to be quasi-long 
range exactly at the lower critical dimension d = 2. In the last case, spectral peaks on the 
oscillating side of the Hopf bifurcation are expected to exhibit power-law tails with non- 
universal exponents. If the formal analogy with the equilibrium critical point found here 
in c? = 4 — e persists in d = 2, we would expect to see features of the Kosterlitz-Thouless 
universality class in systems of coupled oscillators in this dimension. 

Critical oscillators are ideally suited for nonlinear signal detection and amplification. 
Indeed, close to the critical frequency, the linear response function exhibits divergent be- 
haviors, indicative of a high sensitivity of the system. It has been suggested that the ear of 
vertebrates uses critical oscillations of mechanosensitive hair cells for the detection of weak 
sounds and that the properties of the critical point provide the basis to explain the observed 
compressive nonlinear response to mechanical stimuli and to frequency selectivity in the ear 



I, |23[ . The correlation and response functions of single mechanosensory hair bundles have 



been determined experimentally 



2l|. These single cell experiments detected vibrations at 



the scale of tens of nanometers. There, the Hopf bifurcation was concealed by finite-size 
effects but its signature could be observed. In the cochlea of mammals, power-law responses 
over several orders of magnitude have been seen j^. This suggests that in such systems 
a large number of oscillating degrees of freedom operate collectively and bring the system 
closer to true criticality. 

The critical oscillations discussed here can in principle be realized in artificial systems. 
Nanotechnology aims to build functional units on the sub-micrometer scale. Large arrays 
of nano-rotators or oscillators on patterned substrates coupled to their neighbors by elastic 
or viscous effects would provide a 2-dimensional realization of our field theory. This could 
permit in the future experimental studies of the critical phenomena discussed here. 

We thank Edouard Brezin, Erwin Frey and Kay Wiese for useful discussions. 
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FIG. 1: (a) Graphic representation of the propagators Ga/3 and and of the vertex Uap5^5- (b) 
Examples of Feynman diagrams of the theory to one and two loop order. 
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